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ABSTRACT 

fvj ' Context. Silicate carbon stars show the 10 //m silicate emission, despite their carbon-rich photospheres. They are considered to have 

circumbinary or circum-companion disks, which serve as a reservoir of oxygen-rich material shed by mass loss in the past. 
Aims. We present A'-band spectro-interferometric observations of the silicate carbon star BM Gem using MIDI at the Very Large 

^S^ , Telescope Interferometer (VLTI). Our aim is to probe the spatial distribution of oxygen-rich dust with high spatial resolution. 
Qh' Methods. BM Gem was observed with VLTI/MIDI at 44-62 m baselines using the UT2-UT3 and UT3-UT4 baseline configurations. 
Results. The A'-band visibilities observed for BM Gem show a steep decrease from 8 to ~ 10 yum and a gradual increase longward 
of ~10yum, reflecting the optically thin silicate emission feature emanating from sub-micron-sized amorphous silicate grains. The 
differential phases obtained at baselines of ~44-46 m show significant non-zero values ( 70°) in the central part of the silicate 



o 



^ ' emission feature between ~9 and 1 1 /im, revealing a photocenter shift and the asymmetric nature of the silicate emitting region. The 

' observed A'-band visibilities and differential phases can be fairly explained by a simple geometrical model in which the unresolved 

star is surrounded by a ring with azimuthal brightness modulation. The best-fit model is characterized by a broad ring (~70 mas across 

at 10 /jm) with a bright region which is offset from the unresolved star by ~20 mas at a position angle of ~280°. This model can be 

^ ■ interpreted as a system with a circum-companion disk and is consistent with the spectroscopic signatures of an accretion disk around 

fvj ' an unseen companion recently discovered in the violet spectrum of BM Gem. 
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^^ 1. Introduction star drives an outflow from the circum-companion disk, and the 

^^ ' silicate emission emanates from this outflow. 
^— ^ ' "SiUcate carbon stars" are characterized by oxygen-rich circum- 

> ■ steUar material as evident from the prominent silicate emission. Thanks to the recent progress of high-angular resolution 
despite their carbon-rich photospheres (Little-Marenin [19861 observation techniques, there is growing observational evi- 
Willems & de Jong [T9861 . The detection of masers from O- '1^"'^^ for the presence of circumbinary or circum-companion 
H bearing molecules such as OH and H2O from some of the silicate disks around sihcate carbon stars. Our N-hand spectro- 
■ - carbon stars confirms the presence of oxygen-rich circumsteUar interferometric observations of IRAS08002-3803 with the MID- 
gas (e.g., Nakada et al. 1987 ; EngelslgSiT; Szymczak et al.|200l] '"^'"^^'^ Interferometric instrument (MIDI) at ESO's Very Large 
and references therein). Moreover, optical spectroscopic stud- Telescope Interferometer (VLTI) spatially resolved the dusty en- 
ies (e.g., Llovd-Evans lT990l Chan lT993l Ohnaka & Tsuii fT999l l vironment of a silicate carbon star for the first time (Ohnaka 
reveal that silicate carbon stars have I2c/i3c ratios as low as ^t al. |2006), and our radiative transfer modeling suggests the 
-4-5 (classified as "J-type"), which is difficult to explain by the presence of an optically thick circumbinary disk in which amor- 
present steUar evolution theory. Currendy, silicate carbon stars P^ous silicate and a second grain species (amorphous carbon, 
are considered to have an unseen, low-luminosity companion, l^'g^ silicate grains, or metallic iron) coexist. The presence 
The oxygen-rich material which was shed by mass loss when of a circumbinary disk with multiple grain populations is also 
the primary star was an M giant can be stored in a circumbinary suggested for another silicate carbon star IRAS 1 8006-32 1 3 by 
disk (MorrisUMi Lloyd-Evans 1990) or in a circumstellar disk ^eroo et al. ([2003. The spatial distributions of die 22 GHz 
around the companion ("circum-companion disk", Yamamura et H2O maser emission toward two silicate carbon stars, V778 Cyg 
al. 12000). even after the primary star becomes a cai'bon star. In ^n^ EU And, have recently been resolved by Szczerba et al 
flie latter case, Yamamura et al. (I2OOOI 1 suggest that the primary (2006), Engels (priv. comm.), and Ohnaka & Boboltz (200i 

using MERLIN and the Very Long Baseline Array (VLBA). 

The observed spatio-kinematic structures are consistent with the 

Send offprint requests to: K. Ohnaka presence of circum-companion disks. However, the formation 

* Based on observations made with the Very Large Telescope mechanisms of these disks are little understood, and the pos- 

Interferometer of the European Southern Observatory. Program ID: sible link between the disk formation and the anomalously low 

078.D-0292 '^C/'-'C ratios is by no means clear, either. 
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Table 1. MIDI observations of BM Gem: night, time of obser- 
vation (Coordinated Universal Time=UTC), telescope configu- 
ration (Tel.), projected baseline length Bp, position angle of the 
projected baseline on the sky (P.A.), and seeing in the visible 
at the times of the observations of BM Gem and the calibrators 
listed in Table |2] 



# 


Night 


'obs 

(UTC) 


Tel. 


(m) 


P.A. 

n 


Seeing 


1 


2006 Dec. 29 


07:29:47 


UT3-4 


45.7 


93.1 


1'.'0-1'.'5 


2 


2006 Dec. 31 


06:15:26 


UT3-4 


55.2 


99.1 


a'7-i'.'o 


3 


2007 Jan. 04 


06:31:59 


UT2-3 


43.8 


49.8 


ff.'S 


4 


2007 Jan. 05 


04:17:10 


UT3-4 


62.1 


108.0 


O'.'7-l'.'O 



Table 2. MIDI observations of calibrators: 12 /im fluxes (Fn), 
uniform-disk diameters (liuo)- night, and the time stamp (fob.s)- 
The uniform-disk diameters were taken from the CalVin list 
available at ESO (http://www.eso.org/observing/etc/). The data 
set used for spectrophotometric calibration is marked with f. 



Calibrator F12 

(Jy) 



dvD 
(mas) 



night 



Us (UTC) 



HD48433 


8.5 


2.10 ±0.13 


2006 Dec. 29 

2006 Dec. 31 

2007 Jan. 04 


07:52:09 

06:42:34+ 

06:08:26 


HD47205 


9.7 


2.30 ±0.12 


2007 Jan. 05 


03:41:13 


HD50778 


24.6 


3.95 ± 0.22 


2006 Dec. 31 
2006 Dec. 31 


01:35:21 
05:35:01 


HD98430 


19.1 


3.09 ±0.16 


2007 Jan. 04 


07:03:45 



Direct detection of a companion toward a silicate carbon star 
is hampered by the huge luminosity contrast between the pri- 
mary carbon star (lO-'-lO'* Lq) and the companion (^1 Lq), as 
well as small angular separations (e.g., Engels & Leinert 1994). 
However, Izumiura ( 120031 ) and Izumiura et al. (2008 ) have over- 
come the first problem by taking advantage of the fact that the vi- 
olet flux of a carbon star (primary star) is heavily suppressed due 
to molecular and atomic absorption. Using the High Dispersion 
Spectrograph (HDS) on the Subaru Telescope, they discovered 
continuum emission below 4000 A and Balmer lines showing 
P Cygni profiles with an outflow velocity of ^400 km s"' toward 
the silicate carbon star BM Gem. These spectroscopic signatures 
strongly suggest the presence of an accretion disk around an un- 
seen companion, although it could not be spatially resolved, and 
its nature — whether it is a main-sequence star or a white dwarf — 
remains to be determined. Comparison between the observation- 
ally derived accretion luminosity (~0.2 Lq) and the theoretical 
prediction from the Bondi-Hoyle-type accretion resulted in es- 
timated binary separations of ~3-60 AU for a main-sequence 
star companion and ~ 1 00-500 AU for a white dwarf companion. 
These values translate into angular separations of 2.5-50 mas 
and 83-420 mas for a distance of 1.2 kpc adopted for BM Gem 
by Izumiura et al. (2008). 

Given these estimated separations, it is challenging to re- 
solve the circum-companion disk toward BM Gem with a single- 
dish telescope. On the other hand, infrared interferometry pro- 
vides us with a unique opportunity to achieve sufficiently high 
spatial resolution. In this paper, we present the results of mid- 
infrared spectro-interferometric observations of BM Gem with 
VLTI/MIDI to probe the presence of the circum-companion disk 
with dust thermal emission. 



2. Observations 

Our MIDI observations of BM Gem and the calibrators are sum- 
marized in Tables [1] and |2l respectively. BM Gem was observed 
in the HIGH_SENS mode using a prism {A/AA ^ 30 at 10 yum), 
and four data sets were taken using the UT2-UT3 and UT3-UT4 
configurations with projected baseline lengths of 44-62 m. A 
detailed description of the instrument and the observing proce- 
dure are given in Przygodda et al. (120031 ), Leinert et al. ( 120041 ), 
and Chesneau et al. (120051 ). We used the MIA-hEWS package 
ver. 1 .43 to reduce the MIDI data (Leinert et al. 2004; Jaff'e 2004 ). 
While both of MIA and EWS can extract visibility amplitude, 
EWS can also extract the so-called "differential phase", which 



contains information on the wavelength dependence of the ob- 
ject's phase. However, it should be kept in mind that two pieces 
of information on the object's phase are lost in the derivation of 
differential phase as described in Jaffe (2004): the absolute phase 
offset and the phase gradient with respect to wavenumber. 

We assessed the data quality and derived the calibrated visi- 
bilities as described in Ohnaka et al. (I2008I) . Since it is difficult 
to properly estimate the errors of the calibrated visibilities from 
just 1-3 calibrator measurements, we assumed a total relative 
error of 15% as adopted in Ohnaka et al. (2008). The visibili- 
ties obtained with MIA and EWS are in good agreement, and 
we only show the results derived with EWS in the discussion 
below. We also extracted the absolutely calibrated A^-band spec- 
trum of BM Gem from our MIDI data as described in Ohnaka et 
al. (120071 ). A calibrator observed at airmass similar to our target 
was used for the spectrophotometric calibration and is marked 
with t in Table |2] The error of the absolutely calibrated MIDI 
spectrum is difficult to estimate properly with only one calibra- 
tor, so we assumed errors of 10-20% as in Ohnaka et al. (i2007] ). 



3. Results 

Figure [1^ shows the MIDI spectrum of BM Gem, together 
with IRAS Low Resolution Spectrum (LRSJ3 and the photo- 
metric data from the Midcourse Space Experiment (MSX, Egan 
et al. 120031 ) and the IRAS Point Source Catalog (PSC). The 
MIDI spectrum and IRAS/LRS show that the flux level has lit- 
tle changed over the last 23 years, and the flux variation long- 
ward of 13 jum is also absent, as can be seen from the IRAS/LRS 
and MSX data taken in 1996-1997. This stabihty of the mid- 
infrared flux is the same as found for V778 Cyg (Yamamura et 
al. 120001 ). IRAS080 02-380 3 (Ohnaka et al. ^ODS*), IRAS 18006- 
3213 (Deroo et al. 12007] ). and EU And (Ohnaka & Boboltz 
I20081 ). 

In Fig. [TJ), we also plot the spectrum of V778 Cyg obtained 
with the Short Wavelength Spectrometer (SWS) onboard the 
Infrared Space Observatory (ISO) presented in Yamamura et al. 
(12000b as weU as those of IRAS08002-3803 and I RAS18006- 
3213 obtaine d with MIDI by Ohnaka et al. ([2006b and Deroo 
et al. ( 120071 ), respectively. BM Gem and V778 Cyg show a 
well-defined peak at -9.8 jum, while IRAS08002-3803 and 
IRAS 18006-32 13 show broader features. Yamamura et al. 



Available at|http://www.strw.leidenumv.nl/~ nevec/MIDI 



- Downloaded from http://stafF.gemini.edu/"' kvolk/getlrs_plot.html 
with the correction of the absolute calibration by Cohen et al. ( I1992) . 
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Fig. 1. a: A^-band spectrum of BM Gem extracted from the 
MIDI data, plotted together with IRAS/LRS as well as the photo- 
metric data from the MSX and IRAS PSC catalogs. The width of 
each MSX band is shown with the horizontal bar. b: Comparison 
of the A^-band spectra of four silicate carbon stars: BM Gem, 
V778 Cyg, IRAS08002-3803, and IRAS 18006-32 13. c: SED of 
BM Gem constructed from the (spectro)photometric data avail- 
able in the literature (Hipparcos & Tycho Catalogs, ESA I 19971 
2MASS, Cutri et al. |2003t Noguchi et al. 1981 ; COBE/DIRBE, 
Smith et al. 2004; MSX; IRAS/PSC; IRAS/LRS; MIDI). The 
COBE/DIRBE data longward of 12 fim are not used, because of 
the large errors. The dashed line represents the stellar flux con- 
tribution, which is approximated with a blackbody of 3000 K. 



(I2OOOI 1 show that the siUcate feature of V778 Cyg is charac- 
terized by optically thin emission from sub-micron-sized amor- 
phous silicate grains. The close resemblance between BM Gem 
and V778 Cyg indicates the predominance of such small, amor- 
phous silicate grains around BM Gem, and the silicate dust to- 



ward this object is optically thin. This is also consistent with the 
spectral energy distribution (SED) of BM Gem plotted in Fig.[T];, 
which shows only a modest infrared excess longward of ~8 //m. 
As Yamamura et al. (I2000I I argue, such small grains in an opti- 
cally thin environment are blown away by the radiation pressure 
of the primary (carbon-rich) star, suggesting that a circumbinary 
(or circum-primary) disk cannot exist for a long period of time. 
This argument is in line with our interpretation of the A^-band 
visibilities and differential phases presented in Sect|5] 

On the other hand, the broad A^-band spectra of IRAS08002- 
3803 and IRAS 18006-32 13 indicate the presence of grain 
species other than small silicate. In particular, featureless dust 
emission from large grains, amorphous carbon, metallic iron, or 
iron sulfide is suggested for both objects (Ohnaka et al. 20061 
Deroo et al. 120071 1. Also, the SEDs of these two silicate carbon 
stars show much larger infrared excesses longward of 3-4 jum 
(see Fig. 5 in Ohnaka et al. [20061 and Deroo et al. l2007l l than 
BM Gem. This suggests the presence of optically thick disks, 
which is actually confirmed by the radiative transfer modeling 
of Ohnaka et al. (2006). These differences in mineralogy and 
optical thickness seem to affect the wavelength dependence of 
the observed A^-band visibilities and differential phases, as dis- 
cussed below. 

The A^-band visibilities observed toward BM Gem shown in 
Fig. |2^ basically reflect the optically thin silicate emission: the 
visibilities are the lowest at ~10;um, where the flux contribution 
of the extended silicate emission is the highest. The baseline de- 
pendence of the visibilities observed at ^9 yum and at ^1 1 yum is 
not significant in the baseline range from 44 to 62 m, given the 
errors of the measurements. This means that the extended sili- 
cate emission is completely resolved out and that the sampled 
visibility component corresponds to the unresolved central star 
(its angular diameter is estimated to be ~2 mas as explained be- 
low), for which the visibility is constant for the baselines used in 
our MIDI observations. In this case, the observed visibilities rep- 
resent the fractional flux contribution from the star On the other 
hand, the visibilities observed between 9 and 1 1 fim increase 
with baseline length from 44-46 m to 62 m, as can be seen in 
the data #1, #2, and #4 taken at roughly the same position angle. 
Therefore, the observed visibilities contain some information on 
the geometry of the silicate emitting region and will be modeled 
in the next section. 

Figure |2j) shows the calibrated differential phases of 
BM Gem. Significant non-zero differential phases (~ -70° at 
10 ;um) are detected in the central part of the silicate emission 
between ~9 and 1 1 /im in the data sets obtained at the shortest 
baselines (#1 and #3), while the differential phases observed at 
the longer baselines (#2 and #4) are close to zero. As shown in 
Fig. |2};, the differential phases derived from the calibrators for 
the data sets #1 and #3 are close to zero. The phases of two cal- 
ibrators for the data set #3 are in agreement within +5°, which 
we take as the error of the observed differential phases (the phase 
above 12 fim is more uncertain due to the increasing effect of the 
atmospheric water vapor). This confirms that the non-zero dif- 
ferential phases observed toward BM Gem are real, revealing a 
photocenter shift in the silicate feature caused by an asymmet- 
ric circumstellar environment. The non-zero values are observed 
only in the central part of the silicate feature, although one might 
expect the differential phase to deviate from zero over the whole 
feature. However, as explained above, the silicate emission is 
completely resolved out at ^9 fim and ^11 yum even at the short- 
est baselines. The visibilities observed at these wavelengths cor- 
respond to the unresolved, thus symmetric star This is why dif- 
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Fig. 2. a: A^-band visibilities of BM Gem obtained with MIDI. 
The error bars are shown only for the data sets #1 and #4 
for visual clarity, but all error bars are drawn in the inset, 
which shows the enlarged view of the visibilities observed be- 
tween 9.4 and 11.4 fim. The data set numbers are defined in 
Table[T] b: Calibrated differential phases obtained for BM Gem. 
c: Instrumental differential phases as measured on the calibra- 
tors. The open circles represent the phase observed on the same 
night as the data set #1, while the open squares and upside-down 
triangles represent that observed together with #3. 



ferential phases are nearly zero at these wavelengths, as well as 
in the data taken with the longer baselines. 

The wavelength dependence of the A^-band visibilities ob- 
served for BM Gem is in marked contrast to that previously 
observed toward IRAS08002-3803 and IRAS 18006 -32 13 with 
MIDI by Ohnaka et al. (|2006) and Deroo et al. dSOOTl l. re- 
spectively. The A^-band visibilities of these two objects show 



an increase from 8 to ~10 /im and remain nearly constant 
or slightly decrease longward of 10 /urn. We also note that 
Deroo et al. (2007) observed non-zero differential phases to- 
ward IRAS 18006-32 13, which are characterized by a signifi- 
cant phase jump at ~8.3 fim, completely different from that ob- 
served toward BM Gem. Given that the visibility levels observed 
toward IRAS 18006-32 13 range from -0.1 to -0.7, the differ- 
ence in the A^-band visibilities and differential phases cannot 
be solely attributed to a difference in the object's angular size 
and/or baselines used in the observations. As outlined above, 
the dusty environment of BM Gem is very different in miner- 
alogy and optical thickness from that of IRAS08002-3803 and 
IRAS 18006-32 13. Therefore, the distinct A^-band visibilities and 
differential phases observed toward BM Gem most likely re- 
flect the difference in dust chemistry and the structure of the 
circumstellar environment between BM Gem and IRAS08002- 
3803/IRAS18006-3213. 

4. Geometrical model 

To characterize the asymmetric silicate emitting region, we 
adopted an "asymmetric-disk" model, in which an unresolved 
star is surrounded by a ring with a Gaussian radial cross sec- 
tion, and the azimuthal brightness distribution of the ring is si- 
nusoidally modulated. This model was used to interpret the inter- 
ferometric data on Herbig Ae/Be stars by Monnier et al. (:2006l l. 
The intensity at the radial distance r from the star and the posi- 
tion angle (p is expressed as 



I,i(r, (p) oc e 



X {flCOS(0-<l)) + 1) 



where Rx and cr^ denote the radius and the Gaussian width of the 
ring, respectively. The intensity contrast between the brightest 
and faintest points of the ring is determined by the parameter a, 
and O is the position angle of the brightest point. Rx and cta are 
dependent on wavelength, but they cannot be well constrained 
by the current MIDI data at ^9 /im and ^11 //m, because the ex- 
tended component is resolved out at these wavelengths. To take 
into account that the emission at the longer wavelengths origi- 
nates in the cooler, outer region, we assumed the radius to be 
expressed as 



R,^R 



lO/im 



[i^'- 



where Ru)^^ denotes the ring radius at 10 yum. The ratio be- 
tween R;i and cr^ was assumed to be independent of wavelength. 
Therefore, the free parameters of the model are R\Q^m, P, the 
ratio (Ta/Ra, a and O. 

We first derived the ratio of the stellar flux to the total flux at 
each wavelength as follows. If we approximate the stellar radia- 
tion with a blackbody of effective temperature T^ff (3000 K for 
BM Gem, Ohnaka & Tsuji 1999), the observed flux Ff^ (star -i- 
extended dust emission) is given by 

where 0* and /* are the stellar angular radius and the ratio of the 
stellar flux to the total flux, respectively. As mentioned above, 
the visibilities at $9 ;um and ^11 jjm represent the fractional 
flux contribution of the star Using the observed 8 yum visibility 
for /* and the observed 8 jim flux, 6^, was derived to be 1.1 mas. 
Using this 0* and the MIDI spectrum, we determined /* for all 
wavelengths in the A^ band. We note here that as Fig. (Tp shows, 
the SED of BM Gem shortward of ~5 fim is well fitted by the 
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Fig. 3. a: A^-band visibilities of BM Gem obtained with MIDI (symbols) and those predicted by the best-fit asymmetric-disk model 
discussed in Sect. |4] (lines). The error bars are shown only in the inset, which shows the enlarged view of the visibilities observed 
between 9.4 and 11.4 jjm. The data set numbers are defined in Table[T] The model visibilities are shown by the solid, dashed-dotted, 
dashed, and dotted lines for the data sets #1, #2, #3, and #4, respectively, with the same colors as the corresponding data sets, b: 
Calibrated differential phases obtained for BM Gem (symbols) and those predicted by the asymmetric-disk model (lines) shown in 
the same manner as in the panel a. The dashed-dotted and dotted lines (red and orange lines) are almost entirely overlapping, c-e: 
8, 10, and 13 fim images predicted by the best-fit model. The radius of the ring is 20 mas with a Gaussian radial cross section with 
a FWHM of 27 mas at 10 /jm. The azimuthal intensity modulation is characterized by a = 0.2 and <1> = 280°. The image at each 
wavelength is normalized so that the integrated flux is 1 (the color on the central star is saturated). A color version of this figure is 
available in the electronic edition. 



blackbody of 3000 K, suggesting the absence of hot dust, which 
would have manifested itself as an additional (unresolved) com- 
ponent in our MIDI observations. This verifies our assumption 
that the observed visibility at 8 /im represents the fractional flux 
of the star. 

The intensities of the point-source-like star and the extended 
component at each wavelength were set so that the flux ratio 
is consistent with the derived /*. The visibility amplitude and 
phase were calculated from the two-dimensional complex visi- 



bility computed at each wavelength. The differential phase was 
calculated in the same manner as the EWS software derives it 
from MIDI data: a least-square-fit with a straight line (i.e., a 
linear function with respect to wavenumber) was performed on 
the model phase as a function of wavenumber. The differential 
phase was then derived by subtracting the fitted linear compo- 
nent from the model phase. The sign of differential phase mea- 
sured with MIDI is defined as described in Ratzka et al. (,2008] l. 
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In the model fitting, we excluded the observed diff'erential phase 
data above 12 yum because of the larger systematic uncertainties. 
Figures [3^ and ^ show a comparison between the best-fit 
model and the MIDI data. This model is characterized by a broad 
ring with a radius of 20 mas and a Gaussian radial cross section 
with a FWHM of 27 mas at 10 /im (cta/Ra = 0.8). The wave- 
length dependence of the ring radius is characterized by p - 2, 
which results in overall sizes of ~45, 70, and 120 mas at 8, 10, 
and 13 yum, respectively, as shown in Figs. |3};-[3^. The intensity 
of the ring has a mild azimuthal modulation with a = 0.2, and 
the brightest point is offset by 20 mas (at 10 yum) at a position 
angle of 280° with respect to the star Figure [3^ reveals that the 
observed A^-band visibilities are reasonably reproduced by the 
model. As shown in Fig. [3j), the predicted differential phases 
also agree fairly with the observed data. The model predicts the 
non-zero differential phase values to be somewhat too large at 
$9 fim and too small at ^9 yum compared to the observed data. 
Also, the model visibility for the data set #3 is lower than ob- 
served. However, given the simple nature of the model, the fit 
can be regarded as fair The;t'^ value is 3.8 for this best-fit model. 
The formal errors of the ring radius at 10 //m (Rio^m), p, the ratio 
(Ta/Ra, azimuthal brightness modulation (a), and the position an- 
gle of the brightest point (O) are ±2 mas, +0.5, ±0.05, ±0.1, and 
+40°, respectively. However, it should be kept in mind that the 
extended silicate emission is nearly resolved out by our MIDI 
observations, which makes it diflicult to glean information on 
the overall geometry of the silicate emitting region. Therefore, 
the above model parameters would better be regarded as prelim- 
inary despite the rather small formal errors, until they are con- 
firmed by future MIDI observations with a better (m, v) coverage. 

5. Discussion 

Given the spectroscopic signatures of a circum-companion ac- 
cretion disk discovered toward BM Gem by Izumiura ( 120031 ) and 
Izumiura et al. (120081) . it is tempting to interpret the above model 
as a circum-companion disk. In the scenario of Yamamura et al. 
( 1200 0), the radiation pressure from the primary star drives an 
outflow from the circum-companion disk. Depending on the or- 
bital period of the binary system and the outflow velocity, the 
material blown away from the disk may trail like a cometary tail 
and show a spiral structure surrounding the whole binary sys- 
tem. The material surrounding the whole system (corresponding 
to the "cometary tail") can be represented by the broad ring of 
the asymmetric-disk model, while the circum-companion disk 
(the "coma" of the cometary trail) appears as the offset bright 
region. In our model, the offset of the bright feature is depen- 
dent on wavelength as the ring's radius, although the position 
of the companion is independent of wavelength. However, since 
the mid-infrared emission reflects the distribution of dust, not ex- 
actly the companion itself, such wavelength dependence of the 
offset may also be expected. 

If the above model represents the circum-companion disk, 
the offset of the bright feature from the central star, 13-34 mas, 
corresponds to the angular separation between the primary star 
and the companion. This favors a main-sequence star companion 
rather than a white dwarf companion, because the luminosity 
from an accretion disk around a white dwarf at separations of 
13-34 mas (16^1 AU at a distance of 1.2 kpc) is estimated to 
be ~5-18 Lq in the Bondi-Hoyle formulation (Eq. 5 in Izumiura 
et al. 1200 8). This is significantly higher than the observationally 
derived value of ~0.2 Lq (0.03-0.6 Lq, Izumiura et al. 2008). 
However, it should be noted that the Bondi-Hoyle-type accretion 
could overestimate the accretion rate as much as by an order 



of magnitude (e.g., Nagae |2004| and references therein). In this 
case, the predicted accretion luminosity is ~0.5-1.8 Lq, which 
is marginally compatible with the observed range. Therefore, a 
main-sequence star companion seems more likely, although a 
white dwarf companion cannot be entirely excluded. 

6. Conclusion 

Our VLTI/MIDI observations of the silicate carbon star 
BM Gem have revealed that the A^-band visibilities are char- 
acterized by a steep decrease from 8 to ~10 yum and a gradual 
increase long ward of 10 fim, reflecting the optically thin 10 yum 
emission feature emanating from sub-micron-sized amorphous 
silicate grains. Furthermore, we have detected significant non- 
zero differential phases ( — 70°) between ~9 and 1 1 yum, which 
reveals the asymmetric nature of the silicate emitting region. The 
wavelength dependence of the observed A^-band visibilities and 
differential phases is remarkably different from that previously 
observed toward two silicate carbon stars, IRAS08002-3803 and 
IRAS 18006-32 13. This is most likely to reflect the differences 
in mineralogy, optical thickness, and geometry (circumbinary or 
circum-companion) of the disks. 

Our simple geometrical modeling of the observed visibilities 
and differential phases suggests the presence of a bright region 
which is offset from the central star by 13-34 mas (16-41 AU at 
1 .2 kpc) at a position angle of 280±40°. This bright region can be 
interpreted as a circum-companion disk. The derived offset from 
the central star makes a white dwarf companion unlikely (thus 
leaving a main-sequence star companion more favorable), be- 
cause a white dwarf companion at these separations would pro- 
duce too much accretion luminosity compared to that observed. 

However, the lack of detailed information on the geometry 
of the silicate dust distribution makes this circum-companion 
disk interpretation not yet definitive. For example, we can- 
not yet entirely exclude the possibility that the asymmetric- 
disk model represents a circumbinary disk with an asymmetric 
structure. MIDI observations with shorter baselines are indis- 
pensable for putting tighter constraints on the size and shape 
of the silicate emitting region. Also, high-resolution imaging 
at violet/blue wavelengths is important for direct detection of 
the circum-companion accretion disk. Combination of the vi- 
olet/^blue high-resolution imaging and MIDI observations will 
provide definitive confirmation that silicate dust is stored in a 
circum-companion disk. 
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